Renyi entropy of a line in two-dimensional Ising models 



J.-M. Stephan, 1 G. Misguich, 1 and V. Pasquier 1 

1 Institut de Physique Theorique, CEA, IPhT, CNRS, URA 2306, F-91191 Gif-sur-Yvette, France. 

('Dated: June 9th, 2010) 

We consider the two-dimensional (2d) Ising model on a infinitely long cylinder and study the 
probabilities pi to observe a given spin configuration i along a circular section of the cylinder. These 
probabilities also occur as eigenvalues of reduced density matrices in some Rokhsar-Kivelson wave- 
functions. We analyze the subleading constant to the Renyi entropy R n — 1/(1 — n) hi(^]Vp™) and 
discuss its scaling properties at the critical point. Studying three different microscopic realizations, 
we provide numerical evidence that it is universal and behaves in a step-like fashion as a function of 
n, with a discontinuity at the Shannon point n = 1. As a consequence, a field theoretical argument 
based on the replica trick would fail to give the correct value at this point. We nevertheless compute 
it numerically with high precision. Two other values of the Renyi parameter are of special interest: 
n = 1/2 and n = oo are related in a simple way to the AfHeck-Ludwig boundary entropies associated 
to free and fixed boundary conditions respectively. 



I. INTRODUCTION 

The entanglement (or Von Neumann) entropy is in gen- 
eral a difficult quantitypto compute in two-dimensional 
quantum lattice models.tl In Ref. |^ it was however shown 
that for a particular type of wave functions, of type 
dubbed "Rokhsar-Kivelson" (RK), and for particular ge- 
ometries, the calculation simplifies considerably. A lat- 
tice model of statistical mechanics can be used. to define 
a Rokhsar-Kivelson wave-function as follows Ba 
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where the sum runs over the classical configurations and 
E{c) is the energy associated to c (interactions are as- 
sumed to be short-ranged), and the normalization factor 
involves the classical partition function Z. For such a 
state, it has been shown in Ref. || that the eigenvalues 
of the reduced density matrix of a semi-infinite cylin- 
der (with a finite circumference L, see Fig. |l|) are simply 
the classical probabilities pi to observe a given configura- 
tion i at the boundary between A and B. In turn, these 
probabilities can be obtained from the dominant eigen- 
vector of the transfer matrix of the classical model. So, 
the complete entanglement spectrum is encoded in the. 
dominant eigenvector of the classical transfer matrix.Ej 
In this work, we concentrate on the situation where the 
classical model is a two-dimensional Ising model. Each 
probability pi is therefore associated to a given configu- 
ration i of the spins along the "ring" of length L which 
separates the regions A and B (Fig. [I]). Specifically, we 
are interested in the behavior of the Renyi entropies 
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which is the Shannon entropy (or Von Neumann in the 
quantum/RK point of viewd). n |-. 

As discussed in previous studies pu R n {T, L) scales lin- 
early with perimeter L of the cylinder, even at the critical 
temperature. However, the most interesting piece of in- 
formation is the first subleading correction, r n (T). For 
a given temperature T, the later is defined through an 
expansion of R„(T, L) for large L: 



R n (T, L) ~ a n {T)L + r n (T) + o(l) 



(4) 



and is of order one.Ej Contrary to the coefficient a n , r n 
have been argued to be universal. In the case of Ising 
models, r\(T > T c ) = in the high temperature phase, 
and ri(T < T c ) = ln(2) in the low temperature phase.tl 
At the critical point, the previous numerical calculations 
(up to L = 36) lead to n(T = 2V) ~ 0.2544.E3 The nu- 
merical results presented in Sec. O significantly increase 
the precision on this number: ri(T c ) = 0.2543925(5). 
Furthermore, we confirm its universal character by check- 
ing the agreement between three microscopically differ- 
ent realizations of the 2d critical Ising models: on the 
square and triangular lattice, and using the Ising chain 
in transverse field (ICTF). At present, we are not aware 
of any field theory method which is able to compute this 
number. 



In Sec. Ill we analyze the finite-size scaling of ri(/z, L) 
in the vicinity of the critical point, using numerical (but 
exact) calculations for the ICTF. There, the parameter 
(i measures the ratio of the spin-spin interaction over the 
strength of the external magnetic field and plays the role 
of the temperature in the classical Ising model. Away 
from [i = we conclude that ri(/i, L) only depends 
on L(p — 1) in the critical regime, which is consistent 
with a correlation length diverging as l/(fJ, — 1) close to 
the critical point (located at fi c = 1). In particular, we 
confirm the step-like shape of ri(/i,i = oo). 

In Sec. IV we analyze the finite-size scaling of r n (p = 
fx c = 1,L) in the vicinity of n = 1/2 and n = 1, again 
with the ICTF (up to L = 44 sites). The case n = 1/2 
turns out to be exactly solvable (Sec. ^|) and related to 
the "ground-state degeneracy" for a critical Ising model 
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FIG. 1. (Color online) Cylinder geometry with L y 3> L. A 
probability is associated to each spin configuration i of the 
boundary (red circle) between A and B. 



with free boundary conditions, as discussed by Affleck 
and Ludwig.fa In the vicinity of n = 1 the numerical data 
strongly suggests a step-like shape of r n (fi = 1,L = oo) 
as a function of the Renyi parameter n: r„(/i = 1,L = 
oo) = for n < 1 and r„(/i = 1,£ = oo) = In 2 for 
n > 1. This result has some important consequence re- 
garding possible field theory approaches. In particular, 
a singularity at n = 1 would invalidate any attempt to 
compute ri from an analytical continuation to n = 1 of 
the n £ N* result (replica trick). 



II. SHANNON ENTROPY AT THE CRITICAL 
POINT 

A. Square and triangular lattices 

We compute the Shanon entropy R\ using the transfer 
matrix T of the ferromagnetic Ising model. We numer- 
ically diagonalize T (in the full space of dimension 2 L ), 
on the square and on the triangular latticesES for sizes up 
to L = 14 and denote by \L) and \E) the left- and right- 
dominant eigenvectors of T (corresponding to the eigen- 
value with the largest modulus). Then, the probability 
Pi of a configuration i is given by: 
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in the limit of a infinitely long cylinder L y ^> L. 

The results for i?i(T c ), obtained by summing over the 
2 L configurations, are shown in Fig. 0. The linear behav- 
ior, Ri(T c ) ~ L is apparent, as well as the fact that the 
data for the two lattices extrapolate to the same value 
~ 0.254 at L = 0. Although the systems are relatively 
small, it shows that ri(T c ) ~ 0.254 does not depend on 



the microscopic lattice geometry, and is therefore very 
likely to be universal. 

B. Ising chain in transverse field 

As a third microscopic realization of the Ising 2d uni- 
versality class, we study the ICTF: 
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This Hamiltonian proportional to the logarithm of 
the transfer matrix of an anisotropic Ising model on 
the square lattice, with couplings along the y direction 
("time") which are much stronger than in the x direction 
("space") . 

This Hamiltonian is transformed into a free fermion 
problem using the standard Jordan- Wigner transforma- 
tion. The later free fermion problem is then diagonalized 
using a Bogoliubov transformation. The ground-state of 
% is then described as the vacuum of the Bogoliubov 
fermions. The critical point is located at fi = 1. For 
fi > 1 the system is in the ordered phase, with sponta- 
neously broken Z2 symmetry ((a x ) ^ 0), and for fi < 1 
the system is in the disordered (paramagnetic) phase. 

It turns out that the ground-state \G) of the chain is 
simpler to express in a z basis. For an Ising spin config- 
uration \i) labeled by the variables erf = ±1, the proba- 
bility at fi = 1 is : 
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where M is an L x L matrix defined by: 
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This result is derived in Appendix. |A|, where the non- 
critical case fj, ^ 1 is also considered. However, going 
back to the initial 2d classical model, the actual spin 
directions are measured by erf. So, we first compute an 

(z) 

entropy i?„ corresponding to probabilities of z-axis con- 
figurations, and then use the Kramers- Wannier duality 
transforniationta to obtain the desired R n — R^ : 



J?W(M) = ^ ) (l/M)+ln2. 



(9) 



The calculation of R^ amounts to compute 2 L proba- 
bilities, each of which is obtained as a determinant of size 
L x L. Using the translation invariance and the reflec- 
tion symmetry of the chain, the number of,probabilities 
to compute can be reduced to ~ 2 L /(2L)E3 To do so we 
generate one representative for each orbit of spin con- 
figurations (under the action of the lattice symmetries) 
using the "bracelets" enumeration algorithm of Ref. [IJ. 
For the largest size, L = 44 , computing all the proba- 
bilities (2 L = 1.7 x 10 13 ) required about one thousand 
hours of CPU time on a parallel machine. 
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FIG. 2. Shannon entropy Ri(T c ,L) of the Ising models at 
the critical point, plotted as a function of L (a linear term, 
— 0.41L, has been subtracted for clarity). Data for the square 
and triangular lattices and for the Ising chain in transverse 
field are compared. The data are well reproduced by Ri(L) ~ 
aL+ri+b/L and the subleading constant n is evaluated using 
the three largest sizes. Each line represents the leading term 
and the constant, aL + ri. The subleading term appears to 
be the same n ~ 0.254 for the three microscopic models. 



The data for the Shannon entropy R\ are plotted in 
Fig. H and given in Table Q. They significantly extend 
the results published in Ref. pi The columns r[ X \ 

(5) 

and r\ correspond to three different ways to extract 
the subleading constant from R\(fi = with three 

different types of fits (details in the table caption). In all 
cases the result rapidly converges and, using the largest 
size (L = 44 spins) we estimate that r\ = 0.2543925(5) 
at L = oo. 



III. fj, AWAY FROM THE CRITICAL POINT 




FIG. 3. Subleading constant n of the Shannon entropy of the 
Ising chain as a function of fi. The critical point corresponds 
to n = 1. Inset: ri(fi) for different system sizes, plotted as a 
function of (/i — 1)L. 




FIG. 4. 2n Ising models glued together at their boundary 
("Ising book"). In our case, each "page" has periodic boundary 
conditions along the horizontal axis and is semi infinite in the 
vertical direction. Fig. [l] corresponds to 2 pages (n = I). 



In this section, we investigate the behavior of ri in 
the vicinity of the critical point, by considering the Ising 
chain in transverse field away from [i = 1. The results 
are summarized in Fig. ||. 

In this plot, r\[p) is extracted from Ri(L, /i) using a 
fit to ai(fi)L + ri(/i) + b\{p)jL with three consecutive 
values of L. For the size we have studied (here L < 38), 
there is still some visible finite-size effects. In particular, 
the marked oscillations in the vicinity of fi = 1 are not 
converged to the L = oo limit. In fact, it is reasonable 
to expect the curves to gradually approach a step-like 
function as L increases: r\ = for /.t < 1 and T\ = ln(2) 
for fx > 1. 

This scenario, anticipated in Ref. is corroborated by 
the scaling shown in the inset of FigTTft When plotted as 
a function of (/i — 1)L, the data for different system sizes 
and different values of /i collapse onto a single - and very 
likely universal - curve. This can be understood from 
the fact that the correlation length £ of the Ising model 



diverges as l/\fx— 1\ at the transition, and if one assumes 
that n is a function of in the critical region. If 

correct, it immediately implies that ri(n) is a step-like 
function in the thermodynamic limit. 

IV. RENYI ENTROPY AWAY FROM n = 1 

We now consider the effect of changing the Renyi pa- 
rameter n. When 2n is an integer, R n has an interpre- 
tation in terms of the free energy of k = 2n semi-infinite 
Ising models which are "glued" together at their bound- 
ary (see Fig. [|). 

Using the transfer matrix point of view, it is simple 

to see that p^ 2 is (proportional to) the probability to 
observe the spin configuration i on a circle along which 
k Ising models (defined on semi-infinite cylinders) are 
forced to coincide. This was used in Refs. |rj| and |7| in 
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ML) 
7.02789845748593 
7.87432026832476 
8.72076710746883 
9.56723215961776 
10.4137108773778 
11.2602001105626 
12.1066976079502 
12.9532017180203 
13.7997112017585 
14.6462251114521 
15.4927427098430 
16.3392634147881 
17.1857867605076 
18.0323123698967 
18.8788399343835 



0.4232735600 
0.4232735603 
0.4232735604 
0.4232735605 
0.4232735605 
0.4232735606 
0.4232735605 
0.4232735608 
0.4232735600 
0.4232735614 
0.4232735597 
0.4232735610 
0.4232735606 
0.4232735603 
0.4232735602 



' i 

0.2544012149 
0.2543983072 
0.2543965648 
0.2543954570 
0.2543947190 
0.2543942083 
0.2543938437 
0.2543935763 
0.2543933760 
0.2543932227 
0.2543931036 
0.2543930095 
0.2543929343 
0.2543928734 
0.2543928237 



'_i 

0.2543924985 
0.2543925177 
0.2543925180 
0.2543925156 
0.2543925136 
0.2543925110 
0.2543925139 
0.2543925001 
0.2543925306 
0.2543924796 
0.2543925326 
0.2543925037 
0.2543924999 
0.2543925161 
0.2543925300 



0.2543925471 
0.2543925302 
0.2543925183 
0.2543925130 
0.2543925110 
0.2543925072 
0.2543925188 
0.2543924741 
0.2543925939 
0.2543923635 
0.2543926640 
0.2543924262 
0.2543924890 
0.2543925658 
0.2543925757 



TABLE I. (Color online) Shannon entropy Ri(L, fi — 1) of the critical Ising chain in transverse field as a function of the system 



size L. The subleading constant ri is extracted using three different fits: 



using the three following system sizes: L, L — 2, L — 4. 



is obtained by a fit to Ri(L) 



is obtained by a fit to Ri(L) ~ aL + ri + bL + ■ 



: aL + ri + bL 1 
• eL~ A using the 



six system sizes L, L — 2, • ■ • , L — 10. rj is obtained by a fit to Ri (L) ~ aL + ri + bL~ +■ ■ ■ fL using the seven system sizes 
L, L — 2, • • • , L — 12. a is the coefficient of the extensive (and non-universal) term, extracted from the seven-point fit above. 
From this analysis, our best estimate for L = oo is n = 0.2543925(5). 



some field theory calculations, but it is also true at the 
microscopic level. The interpretation above does not ap- 
ply when 2n is not a positive integer, but R n (L,fi) can 
still be computed numerically for any n > 0. 



with ours for n > 1, but not at n = 1. 



B. 



A. Renyi parameter n = 2 and above 

When n goes to infinity, only the spin configuration 
with the largest probability contributes to R n . For the 
ferromagnetic Ising models we consider (including the 
quantum chain in transverse field), this configuration is 
two-fold degenerate and corresponds to a fully polarized 
ferromagnetic state, | tt ' ' ' t) or I -14 ' ' ' D ■ I 11 other 
words, taking the limit n — > oo amounts to study a semi- 
infinite Ising model with ferromagnetic boundary condi- 
tions. The corresponding probability, p ma x, [behaves as 
— ln(p max ) ~ aL + ln(2) at the critical point.tl The sub- 
leading constant, ln(2), is nothing but (twice) the u g fac- 
tor" associated to this conform ally i nvariant boundary 
condition (more details in Sec. 1VB ). This implies for 
the Renyi entropies that the subleading constant r n (fi c ) 
is ln(2) at n = oo. 

In fact, for n > 2 the Fig. |^ shows that even relatively 
small systems give r n ([i c ) very close to In 2. Table [fj is 
an analysis showing that r2(/i = 1) ~ In 2 with a great 
accuracy, of the order of 10 -8 . Since the convergence to 
In 2 is even faster when n > 2, there is practically no 
doubt that r„(/^ c ) is exactly In 2 for n > 2.E3 

As a consequence, an analytical continuation of this 
result to n = 1 would erroneously give ri(fi c ) = In 2 (in- 
stead of 0.25439). In particular, we note that the results 
of Ref. ^ (which use a replica technique) are in agreement 



The special value n = ^ corresponds to the free energy 
of a single Ising model defined on a semi-infinite cylin- 
der (keeping only part A in Fig. |]), and can be treated 
exactly. Using the transfer matrix point of view, it is in- 
deed simple to see that v /p7 is proportional to the prob- 
ability to observe the spin configuration i at the edge of 
a semi-infinite Ising model (contrary to pi which is the 
probability to observe i in the bulk). 

As far as the universal properties are concerned, we 
can study the ground state of the quantum Ising chain 
(Eq. (^|)) rather than the transfer matrix of the classical 
2d model. Denoting by \G) the ground-state of the chain, 
we have ^Jpl = (G\i) and the Renyi entropy Ri can be 
written as: 



R i (L)=2]n[{G\ J2 l°f-°i>l (10) 
= 21n(2 L / 2 (G|frcc)) (11) 



where |free) = |{of = 1}) is the state where all the spins 
point in the z direction. It turns out that the latter state 
is the vacuum of Jordan- Wigner Fermion and that the 
scalar product in Eq. ( [Tl] ) can be obtained as a particular 
case of Eqs. (^0). At the critical point (/x = 1), the result 
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FIG. 5. (color online) Subleading constant r n of the Renyi entropy of the Ising chain (ICTF, at [i = 1) and classical Ising 
models (triangular and square lattice, at T = T c ) as a function of the Renyi parameter n. The (slow) convergence towards 
a step function can be observed. Inset : when plotted as a function of (n — 1)L ' 25 , the data collapse reasonably well onto a 
single curve. For each value of n, r n is obtained by fitting the data for R„(L) to ~ aL + r n + bL^ 1 + cL~ 2 using four system 
size : L, L — 2, L — 4 and L — 6, as indicated. 
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20 
28 
36 
44 


4.95205232373074 
6.66741530818944 
8.38060934985332 
10.0928119559937 


0.2138075040 
0.2138074244 
0.2138074200 
0.2138074203 


0.9989748222 0.9999928713 0.9999877126 
0.9996525352 0.9999971726 0.9999989449 
0.9998432968 0.9999991645 0.9999998996 
0.9999165184 0.9999996643 0.9999997718 



TABLE II. Renyi entropy R2(L,fi — 1) of the critical Ising chain in transverse field as a function of the system size L. 
The subleading constant r2(fJ, c ) is extracted using three different fits (same as in Table Q): ri, 1 ' is obtained by a fit to R2(L) ~ 
aL + ri + bL^ 1 using the three following system sizes: L,L — 2,L — 4. is obtained by a fit to R2 (L) ~ aL + ri + bL^ 1 + ■ ■ ■ eL~ 4 
using the six system sizes L, L — 2, • • • , L — 10. r 2 is obtained by a fit to R2(L) ~ aL + n + bL^ 1 + • ■ • fL~ 5 using the seven 
system sizes L, L — 2, • • • , L — 12. a is the coefficient of the extensive (and non-universal) term, extracted from the seven-point 
fit above. These data show that n converges to In 2 (limit L — > 00). Similar results, with an even faster convergence, show that 
r n >2(^c) = ln(2). With the present systems sizes and the present machine accuracy, adding more terms in the 1/L expansion 
does not increase the accuracy on n. 



is particularly simple!! Renyi entropy: 

1 L/1 1 

(G[free)= TT cos (2j + 1)7r , (12) Rl (L, y, = 1) = L In 2 + 2 V lncos ^l±3l. (13) 

3=0 AL 2 ,=0 4L 



and leads to the following exact expression of the n = -|- 



Finally, an Euler-Maclaurin expansion gives the desired 
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finite-size scaling, with a vanishing constant ri 



1) = aiL + ri + o(l) 

' 2 2 v ' 

2K 
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(14) 

(15) 
(16) 



where K ~ 0.91596559 is Catalan's constant. 

The constant term in — In (G| free) has already been 
studied in Ref. The situation where \G) is the ground- 
state of an antiferxQmagnetic spin-^ XXZ chain has also 
been consideredJjO Such a scalar product is closely re- 
lated to the notion of quantum fidelityO In terms of 
a classical 2d Ising model, — 7Tn(G|frcc) is the bound- 
ary contribution to the free energy of a semi-infinite 
Ising model with free boundary conditions imposed at 
the edge. At the critical point, this is a well understood 
quantity from boundary CFT, and the subleading con- 
stant ri corresponds to — 21n<7, where g is the "ground- 
state degeneracy" discussed by Affleck and Ludwig.i In 
the present case of tJae Ising model, n = is in agree- 
ment with <7f rec = l.tffl This result has also been checked 
numerically in Ref. B2J 



C. Critical behavior in the vicinity of n — 1 

The results concerning the subleading constant r n (p = 
1) are summarized in Fig. ||. The behavior of r„(/x = 1) 
has some similarity with that of ri(fi): the curves inter- 
polates between and In 2 with a slope at re — 1 (resp. 
fi = 1) which increases as a function of the system size. 
Here again, it appears that the data for different values 
of n and L collapse onto a single curve when plotted as 
a function of (n — l)i 25 (inset of Fig. ||). The error bar 
on the exponent 0.25 are unfortunately large and difficult 
to estimate, but it indicates (a rather slow) divergence of 
the slope dr n /dn\ n= i when L increases. r„ has also been 
computed for the classical Ising mo dels on the square 
and triangular lattices, as in Sec. II A. The inset of Fig. || 
shows that the r n obtained from the corresponding trans- 
fer matrix calculations are in good agreement with those 
calculated from the ground state of the ICTF. This is a 
strong indication that, in a scaling region around n = 1, 
r n defines a universal curve. The analogy between the 
effects of /i and n suggests that n—1 is a "relevant" pertur- 
bation: going slightly below (resp. above) n = 1 induces 
a drastic change in r n (fx = 1), which immediately (when 
L = oo) goes to (resp. In 2), as in high (resp. low) 
temperature phase of the 2d Ising model. 



D. Vicinity of n — | 

The value n = \ can be treated exactly, as explained 
in Sec. |VB|. However, the free-fermion calculation does 
not extend away from re — i. Still, at n — 0.5 we observe 
(numerically) a crossing of the curves corresponding to 
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FIG. 6. (color online) Subleading constant r n of the Renyi 
entropy of the Ising chain (at the critical point /i = 1), in the 
vicinity of n — 0.5. Inset : subleading constant multiplied by 
L 6 as a function of (n — 0.5). 



different values of L (see Fig. ||). This phenomenon, also 
observed at re = 1, is reminiscent of a critical behavior, 
where the deviation away from re = 1/2 would play the 
role of an irrelevant perturbation away from a fixed point. 
The data can also be collapsed onto a single curve, when 
the y-axis is multiplied by a factor w L 6 . However, 
contrary to the case n—1, this result indicates a reason- 
ably fast convergence towards r n = in the vicinity of 
re = 1/2 (see the inset of Fig. 0). 



V. DISCUSSION AND CONCLUSIONS 

In the present Ising models, r n {fj) seems to take only 
three discrete values. For example in the critical case, we 
find : 



r n ([i = 1) 





0.2543925(5) 
In 2 



n < 1 
n = 1 
n > 1 



(17) 



This is quite different from other models described in 
terms of a free field compactified (with radius R) in the 
long-distance limit. In that case, which is better under- 
stood from a field theory point of view, the system de- 
scribes a line of fixed points and the subleading constant 
r n (R) continuously varies along that critical lineE3 



r n (R) = lnii- 



ln j 



2(n-l) 



(18) 



We have discussed how the special values n = g, and 
re = oo are related to the g factors associated to free and 
fixed boundary conditions of the Ising model. But so far, 
we do not know how to understand r n —i for the critical 
Ising model using CFT. This is certainly an interesting 
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question for future studies. Th is is all th e more chal- 



lenging, as it seems (Sees. 1VA and [V C ) that replica 
methods for the Renyi parameter n are not applicable to 
the Ising critical point for this quantity. 

It is tempting to conjecture that crossings for r n (L) are 
observed whenever the underlying probabilities, ~ Pi{^) n 
describe a conformally invariant setup. It is indeed the 
case at n = 1/2 (Ising boundary with free boundary con- 
ditions), but it is also realized for n = 1, since it corre- 
spond to the bulk probabilities. It would be interesting 
to check this idea on other models, and to investigate 
the likely connection with thfi-theory of line defects in 
conformally invariant systems .EjEil. 
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Appendix A: Probability of a spin configuration 



We consider an Ising chain in transverse field 



L-l 



L-l 



(Al) 



3=0 



3=0 



We assume L to be even, as well as periodic boundary 
conditions u x L = crfi. We wish to find the ground state 
\G) of this Hamiltonian, and to compute all of his com- 
ponents in the basis of the eigenstates of the <rj . 



1. Diagonalization 

As is well known, H can be expressed in terms of free 
fermions, using the Jordan- Wigner transformation : 

&j + ia? = 2c] exp \vk 4 c lj (A2) 



a] = 2c)c 3 - 1, 



(A3) 



where the c, satisfy the canonical anticommutation re- 
lations {cj,c\} = 5j£. This allows to write the Hamilto- 
nian as a quadratic form 

L-l L-l 

n = -^& j c j -l)-^(c]-c j )(c* j+1 +c j+1 ), (A4) 

3=0 3=0 

where the fermions are subject to the following boundary 
condition : 

L-l 

c\ = - exp (ittjV) Cq , Af = c[cj. (A5) 



The parity operator V commutes with % 

L-l 
3=0 



(A6) 



and because of Perron-Frobenius theorem, the ground- 
state lies in the sector V = +1. Therefore, fermions are 
subjected to antiperiodic boundary conditions c L = — cj. 
H can finally be diagonalized by a Bogoliubov transfor- 
mation : 



cos 6id 



isin 9kd}_ k 



(A7) 



k e {{21 + \)tt/L \-L/2<l<L/2- l} (A8) 
fi sin k 



sin 29k = 
cos 26^ = 



a/1 + 2/i cos k + /i 2 

1 + (j, cos k 
yl + 2//cosfc + yu 2 



(A9) 
(A10) 



The new fermions operators dk , d\ satisfy the necessary 
anticommutation relations, and diagonalize H ■ 



Ek = 2\J\ + 2/icos/c 



(All) 
(A12) 



Ek > ensures that the ground-state \G) is the vacuum 
|0) of the d k . 



2. Exact formulae for the spin probabilities 

We define PJ as the projector onto the |cr=±l)^ state: 

Pt = c\ C 3 , r,r[. (A13) 

Pi is then given by 

Pl = p(a , . . . , = (0\PSP? • ■ ■ P L CT -i|0>- (A14) 

Using Wick's theorem, this correlator reduces to a Pfaf- 
fian. To compute it, we need to calculate the four types 



of contractions (c|q), (cjc\), (eje]), (cjC(), which can be 
done using Eq. dAj]). It is worth noticing that all these 
correlators are real in this particular model. We write a 
generic projector as : 

Pj = f]fh (A15) 
with /! = Cj for a = +1 and /j = Cj for a = — 1. Then : 

P? = </oVo/J/i /Li/l-i) 2 (A16) 

= </o t /i t ---/l_ 1 /o/i--./L-i> 2 (A17) 
f2 / A B 



= PF 



-B -A 



(A18) 



8 



where A is antisymmetric, B is symmetric, and Pf de- 
notes the Pfafhan. The matrix elements of A and B arc 



(A19) 
(A20) 



B J( = {f}h) 

Using the relation Pf 2 = det, Eq. (A.18) simplifies into 

A B 



Pi = det 
= det 



—B -A 
A + B B 







B — A 



(A21) 
(A22) 



Eq. (A22) follows from Eq. (A21) by adding the second 
column to the first, and then the first row to the second. 



Finally, 

Pi = dct(A + B) = det M, (A23) 
where M is a L x L matrix with elements 

M jt = (f}(fl+m (A24) 

= \ S » + tzT, cos [ fc (i -Q + 20k] (A25) 

k 

At the critical point (// = 1), 6k = k/4 and the matrix 
elements simplify even further : 



Mjt - 2 Sjt + 2Lsm[7r(j-e+l/2)/L] ' 



(A26) 
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This quantum/classical correspondence works in a rather 
straightfoward way for simple constrained models (such as 
dimer models or vertex models). For other models, such as 
the Ising model considered in this paper, some additional 
care is needed to define the geometry of the A/B bound- 
ary at the microscopic level. In the particular case of 2d 
classical Ising models, the spins living at the frontier be- 
tween A and B have to be "duplicated" to insure that the 
decomposition induced by the classical spin configurations 
is indeestia proper Schmidt decompostion of the RK state. 
See RefB for more details. 

In the quantum point of view, where one studies the en- 
tanglement in a RK wave-function, the dominant (~ L) 
contribution is the boundary (alsorpalled "area") law. 
The field-theory prediction of Ref.Q is n(T = T c ) = ln(2) 
and does not agree with our numerical calculations n Re- 
mark: n corresponds to Sq^' in the notatip«p-pf Ref.B. 
T c = 2/ln(l + V2) on the square latticdl3B and T c = 
4/ln(3) on the triangular latticed 

We also use the property that, for periodic boundary con- 
ditions, n, of = 1 in the ground state \G). Since we work 
in the a z basis, this reduces by another factor two the 
number of probabilities to comp ute. 

In fact, the analysis of Sec. IV C suggests that rn(/x c ) flows 
to its ferromagnetic boundary condition limit, ln(2), as 
soon as n > 1. 

See Eq. (77) in Ref. | The following normalization of the 
radius is used: R = 1 for free fermions, and R — spl at 
SU(2) symmetric (self-dual) point of the XXZ chain. This 
corresponds to a lagrangian £ = i (d^ip) 2 , where the field 
is compactified on a circle of radius R : ip — if + 2nR. 



